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Abstract 



Analytical solutions of precise equations that describe the rf-coupling of two 
cavities through a co-axial cylindrical hole are given for various limited cases. For 
their derivation we have used the method of solution of an infinite set of linear 
algebraic equations, based on its transformation into dual integral equations. 

1 Introduction 

In the papers @, ^, H, we derived precise equations, describing the rf-coupling of two 
cavities through a centerhole of arbitrary dimensions. On the base of these equations we 
numerically calculated the relationship of coupling coefficients versus different parameters 
(frequency, hole radius, etc.). This paper presents analytical solutions of these equations 
for various limited cases. In particular, it is explicitly shown that in the case of small holes 
(a — > 0) the formulated equations agree with those derived in the papers [|, [5], on the 
base of quasi-static approach. Besides, expressions are derived for coupling coefficients 
which are valid up to the second order in the relation of the hole dimension (a) with the 
free-space wavelength (A). For derivation of these expressions we have used the method 
of solving of an infinite set of linear algebraic equations, based on its transformation into 
dual integral equations. 

2 Problem definition. Original equations 

Let us consider the coupling of two cavities through a circular hole with the radius a in a 
separating wall that has the thickness t. For simplicity's sake, we will consider the case of 
two identical cavities, with h -being the cavity radii and d — their length. In the papers 
Jj], H, [| it was demonstrated that if the field is expanded with the short-circuit resonant 
cavity modes and S io- m °des are selected as fundamental, the precise set of equations 
will consist of two equations for the amplitudes of E 010 -modes, where coupling coefficients 
are denned by the way of solution of an infinite set of linear algebraic equations. Let us 
generalize the case considered in |l], Q [|, choosing as fundamental E Qqp -modes of closed 
cavities (q is the number of field variations across the radius, p is the number of field 
variations along the longitudinal coordinates). Using the method, similar to the one in 
[U, @, [| , one can show that the set of equations, describing the system under consideration, 
has the form; 






where 



ry 2 11 
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a q 



^jp is the amplitude of E 0qp -mode in the i-th cavity (i = 1,2). The normalized coupling 
coefficients A$ are determined by the expression: 

oo 

a, = Hu) = j 2 (p q ) E u>?V (a 2 - e 2 ) , (3) 

s=l 

where wy) are the solution of the following set of linear equations: 

+ E 0«/£ } + ^ 2) /i 2) ) = 3tt/«/ (A^ - fi 2 ) , (4) 

s 

^ + E G«m H 2) /i 1} + ^/i?) = Stt/W/ (A 2 ; - fi 2 ) , (5) 



fU) _ J c %m) - ch(q m t/l), j = 1 

Jm s%m) I ch(2q m d/l) - 1, j = 2 



J 

q m = Hrnl/a, l — 2d + t, fi m = ^\ 2 n — fi 2 , = cua/c, fi 2 = fi 2 — n 2 a 2 p 2 /d 2 . 

G m ,s B m s — § m s cth( /irrt) "I - 

27ra 2 ^ 3 J 2 (^) 



(6) 
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d& e pX9 (A^ - £ 2 )(A 2 - 0^ + ^oa^/rf 
a ^ 6 2 Jl(6f)R£ 

bhxi{XL-e}){\ 2 -e}) 



Re 



6 e = X e a/b, Xi = rr\ e J 2 (X e )/2, u e = y/flf - ft 2 , 

9 q cth(v q d/a)/v q - 2ad q / \e p d [v 2 + 7t 2 a 2 p 2 /d 2 ) } , £ = q, 
Oecth(vid/a)/ve, £ ^ q. 



The coefficients have a simple physical sense. Really it is easy to show that the 
tangential electric field component in the left cross-section of the coupling hole E)~\r) 
has the form 

4~\r) = E^r) - ES(r) = E^ p Q^(r) - ^ p Q (2) (r), (9) 

where 

Eop q is the value of the longitudinal (perpendicular to the hole) electric field of (0,q,p)- 

~ (2) 

mode in the first cavity on the left coupling hole cross-section at r = a, while Eq P(1 is the 
same value for the right-hand cavity on the right coupling hole cross-section at the same 
radius. 

From the expression (|9]) it follows that the tangential electric field component on the 
left coupling hole cross-section[| is equal to the difference of two induced fields, each of 

2 The same is true for the right cross section 
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which is proportional to the perpendicular electric field components of Eo, q , p -modes, taken 
to be fundamental. There, the coefficients are the ones of expansion of the appropriate 
functions with the complete set of functions {Ji (X s r/a)}. 

Note that the coefficients w$ can be re-defined which will cause changes in Eqs.((|,|5]), 
in the form of relationship (^), and, consequently, in the above nature of proportionality. 
For example, in @ one can obtain proportionality E^ d to the longitudinal electric field 
component of (0, q, p)-mode at r = 0. While defining w$, we proceeded from the following 
condition: for the two different cavities, in the case t = (infinitely thin wall) and taking 
£"0,5,0-1110 des to be fundamental: Q^(r) = Q ( r ) = Q( r )i then we have in this case 
E!f~\r) = E^ + \r) = Q(r) (^E^^ — E^ ^ . This condition determines the tangential 
electric field component on the hole while having Exi^p-modes as fundamental. With such 
a normalization of w W , their determining set of infinite linear equations (|4],|5|) acquires the 
most symmetrical form. 

Thus, the two-cavity coupling problem, rigorously formulated on the base of the elec- 
tric field expansion with the short-circuit resonant cavity mode, is reduced to the induced 
field definition on the right and left cylindrical hole cross-section. 



3 Infinitely thin wall case 

An important role in the problem of cavity coupling plays the case of infinitely thin wall, 
dividing the cavities (t = 0). In this case, from Eqs. ([|,|5|) it follows that w$ = w!ff = w m . 
In this case the set of equations for w m will take on the form[]: 

£ w s B m , s = 3tt/ {2 (a^ - nj) } . (10) 

For the case t = Ai = A2 = A, where 

A = J 2 o (O q )E^/{x 2 s -0 2 q ). (11) 

3.1 Small coupling hole case (a — ► 0) 

If in Eqs. ( |r0"| , |TTD the hole radius tends to zeroQ then Eq.fllTiD will become: 

s ^ s Jo (xi - e*) (\i - 0*) 2X1' 1 ] 

In order to get the solution for Eq. (|T2"D we will introduce an integer odd function fi(z) 
the values of which in the points z = X s are equal 

f 1 (X s )=w s J 1 (X s ). (13) 



We have neglected terms of order a 5 in the expression (^) for G m ,s 
In this case, as follows from Eqs.(|l],|2|), the coupling coefficients will be proportional to 
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Let us assume that at | z |— > oo /i(z) grows not faster than exp(z), then, in accordance 
with Cauchy theorem, the function (f(z)/Jo(z)) can be expanded in a patial fraction 
series 

oo 

h{z)/Mz)=2zY.Wn/(\l-z 2 ). (14) 

71=1 

Using (PI) , and, also, multiplying Eq.(|T2]) by Ji(X m x)/ Ji(X m ), where < x < 1, and 
doing summation over sub-index m, we will get 

/ f 1 (z)J 1 (xz)dz = &Kx/2, 0<x< 1. (15) 

JO 

By multiplying ([14]) by zJ\(xz) and integrating over 2; from to 00, we will obtain (see 
the Appendix) at x > 1: 

/■oo 

Lx(x) = / zf 1 (z)J 1 (xz)dz = 0, x > 1. (16) 
Jo 

In this way, the set of linear algebraic equations (0) with a complicated coefficients 
matrix that cannot be expressed via elementary functions and can be calculated only 
numerically, has been reduced to two integral equations flI5]),(|T6]). Having determined the 
kind of function fi(z), there is no need in calculating the sum (]TT[) , since 



A = 5>./A2 = Mm AC*)/ (2*J„(*))- (17) 



The method of solving the dual integral equations of the type ( |I7||T6D on the base of 
the Mellin transformation, as well as the property of Cauchy-type integrals, can be found 
in 0. The brief summery of their solutions is given in ||. We shall dwell briefly on a 
simpler method of resolving the system, because it will be used in Sec. 3 for the analysis 
of the infinitely thick wall case. 

Since fi(z) is the odd function it can be represented in the form 



h( z ) = / sm(zt)ri(t)dt. 
Jo 

Substituting this expression in Eq. flTB]) we obtain such integral equation for rj(t): 

rj(t)dt 



, 0, x > 1. 

Vt 2 - x 2 

The solution of this equation is r](t) = for t > 1. Consequently, any function of the type 

fj z ) = f 1 sm( zt)r}(t)dt (19) 
Jo 

satisfies Eq.flTEp. Substituting (|I9"D into ([15]), we obtain the first kind Volterra equation 
Abelian type 



trj(t)dt 37r 2 



y/x 2 - t 2 2 



x, < x < 1, (20) 
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the solution for which can be found in the analytical form. Omitting the intermediate 
formulae, we shall give the final expression for the function fi(z) 



h{z) = 1 {sin(z) - * cost*)} « 2z(l - Q (21) 

The normalized coupling coefficients, as follows from (|17D, is equal to A = 1. Since 
w s = /(A s )/Ji(A s ), then, from (|9p, we will obtain 

El-Hr) = E ^ {r = 0) ~ ^ = Q) . r (22) 
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Thus, on the base of a rigorous electrodynamic description of the two cavity coupling 
system we are the first to prove, by the way of the limit transition a — > 0, the correct- 
ness of the equations formulated in the papers |4], |5], |6| on the basis of the quasi-static 
approximation, and to obtain the expression for the tangential electric field on the hole. 

3.2 The case of small, though finite, values of coupling hole 
radius 

The above method presents the opportunity to obtain analytical expressions for the nor- 
malized coupling coefficients with an accuracy on the order of (a/A) 2 . If a/A is small, 
though finite, then, the coefficients w s in (|TDJ) will be dependent on the hole radius value 
a: w s = w s (a). Let's introduce the function of two variables: 

^(a,z)=2zJ (z)±^- 2 (23) 

n=l A n - z 

We will assume that relative to the variable z the function ip(a, z) will obey the conditions 
formulated in Subsec.2.1. Using the technique, similar to that described in Subsec.2.1, 
the set ( |TD| ) can be reduced to: 

g W) ^ ^ g e) = Q ^ 1<x< fe/0) (24) 
e=i X£ 

a ~ J l (x6 e ) 'ZixJ 1 (xO*) 

^tL ip{a,8 e )R e = — — , < x < 1. (25) 

o £ =1 Xi is* Jo I"* ) 



Letting a — > in Eqs.(|2^),(p5|), we derive a set of equations (|15| , [16D , and, consequently, 
^(0, z) = fi{z), where fi{z) is determined by Eq. fl2"T|) . Let's represent ip(a, z) in the form 

^(a, z) = ip(0, z) + a 2 ip(a, z), (26) 

where (p(a, z) is a function which has the same conditions imposed upon that ip(a, z) does. 
From (|2^ ) , (|25|) it follows that <p(Q,z) satisfies the following equations 



QJ X (x6) (p (0, 9) d6 = 0, x > 1, (27) 

rod 

/ J 1 (x6)(p(0,6)d6 = F(x), 0<x< 1, (28) 
Jo 



where 



Fix) = lim — 

a^O a 2 



The coefficients Ri can be represented as: 



=1 



ft 2 . 
2(91 

It can be shown that the following estimations are true 

a ^ Ji (x9f) , . n , 3n „ / ~\ 
7T- £ ^-^^ (0, fc) = -yx + O a 3 



Then 



= — x 



n 2 - n 2 2Q 2 - n 2 



Aa 2 



-x 



The solution of Eqs.fl27.28) has the form 

n 2 - n 2 



2tt 2 - Q 2 



Aa 2 J1V ' 2a 2 
where f\(z) is determined by the formula (|21|) , while /^(z) is 

, . (3z 2 — 6) sin(z) — z (z 2 — 6) cos(,z) z 
h{ z ) = -7 ~ 7- 



(29) 



The function ip(a, z) ~ ?/>(0, z) + a 2 <^(0, 2;), accurate to the order (a/A) 2 , has the form 



ip(a, z)»ll+ I fi(z) f 2 (z). 

The normalized coupling coefficients A is determined by the relationship 



(30) 



A = Jn 



; A «)4^)/( 2 ^) 



1 (a 



3 (a 



1 fa n 2 
5 V6 "V 20 Vc^'V "20 U W 



1 - 3 ( ~-r\„\ — — I —W. 



(31) 



For the case uj ~ w 9>p the expression (|31|) agrees with that for the generalized polarizability, 
obtained in || at b — > 00 via the variation technique. Note that the expression (fH]) is true 
for the frequency uj that is not close to the resonant frequencies of the non-fundamental 
modes of closed cavities: uj 7^ uj m ,n, if (m,n) 7^ (q,p). 

Knowing ip(a,z), and, consequently, w s (a) = ip(a,X s ) / J\ (A s ), the form of the tan- 
gential electric field around the hole can be reconstructed 



E (- )( , = gSU = 0) - E o,U r = 0) 



(2) 



7T 



X 



1 - 



4 \b 



A + 



Q 2 - 2Q 2 



12 



+ 



2ft? - Q 2 r 



6 



1- 



(32) 



Since in our approach the case Q 
calculation, then, the formulae (pTT 



corresponds to the quasi-static method of field 
) at O = present the solution for the appropriate 
quasi-static problem up to the second-order approximation in (a/b) in the case when the 
"far" field (see, f|,f|) is not homogeneous. 



4 Infinitely thick wall case 

The above analytical results pertain to the case of infinitely thin wall, when the singularity 
at the rim of the hole has the form (r — a) -1 / 2 . It is of interest to consider the possibility 
to apply the above method to a non-zero thickness wall, when the singularity at the rim 
of the hole has the form (r — a) -1 / 3 fl~0 |. The simplest problem in this class, although 
representing an important application, is research into the coupling of a cylindrical cavity 
with a co-axial cylindrical waveguide of the radius a < b. Such a system can be studied 
upon the base the above equations in the limit case t — ► oo: 

4 a 3 

e p Z q>p a q>p = -^ 9 ____ ASip) (33) 
where the normalized coupling coefficients at a — > is equal to 



A = J>,/A 2 , (34) 



while w s are the equation solutions 



2X m + ^ 1 Ws J (\ 2 -e 2 )(xi-9 2 )- x 2 m {6b} 



The set (^5[) is different from the above-studied (IT4 ) in additional addends in the diagonal 
matrix elements. Introducing a function of the type (fT4]), we obtain the following set of 
equations: 

poo 

zf(z)J 1 {xz)dz = 0, x > 1. (36) 



^ fjXm) Ji (xX m ) 1 f°° f( w ( v , 37ra; 

2^ — \ WTV~\ + o / fi. z )Ji{xz)dz = — -, < x < 1. (37) 



m=0 



A m <^l (^m) 2 Jo 



As indicated above, from (|3q) it follows that f(z) should be sought in the form of (|19|) , 
then, from (|37D we get the fact that rj(t) has to be the solution of the Fredholm equation 
of the second kind 

f \ 1 4 f 1 / ,x ^ sin(A m t) sin(A m u) 

vw + - / v(t) }^ — 7 fin \ — dt = 12u - ( 38 > 

7T JO m=1 X m J 1 (X m ) 



8 



Since the kernel of this integral equation is degenerate, then its solution have the appear 
ance 

sin(A m u) 



4 00 
r)(u) = 12m 



where 



vr m=1 X m Ji (A 



rj(t) sm(X m t)dt 



(39) 



the constant coefficients which are the solution of the infinite linear set of equations to be 
easily obtained by way of the appropriate integration Since c m = f (A m ) = w m J\ (A n ) 
(see (p~9|) ) , this system can be represented as: 



w m + - 2^ 



7T 



^ A n Ji (A n ) Ji (A m ) 



sin(A n - A m ) sin(A n + A r 



An — A r 



An + A r 



12 {sin(A m ) - A m cos(A m )} 
A^Ji (A m ) 



(40) 



Thus, for fields with the singularity at the rim of the hole of the type 



- a)- 1 / 3 

we have obtained an analogous system (|40| ) instead of the initial linear set of equations 
(pop. Comparing ( fiTf ) and fl35|), we can see that in the system (|35| ) the matrix coefficients 
are expressed through the improper integrals, whereas in the modified set ([40]) they are 
determined by well studied functions. This considerably facilitates both analytical studies 
and numerical simulation. 

We have carried out numerical calculations of the normalized coupling coefficient A 
determined by the formula fl34[) both on the base of (|35|) and (|40|). For a 200*200 matrix, 
the calculations based on (|35| ) gives A = 0.85835, based on (fiTf ) — A = 0.85854, that 
is the results agree up to 2 x 10 -4 . This result confirms the correctness of the above 
analytical method. Note that the obtained A- value corresponds to the results of a purely 
static analysis [|IT[] carried out on the base of the variation technique for the infinite plane 
with a cylindrical hole. 



5 Conclusion 

On the base of our method of reduction of the infinite linear algebraic equation set to dual 
integral equations, we obtained, in different limited cases, the rigorous analytical solutions 
regarding the two-cavity coupling problem. Alongside with general theory significance, the 
obtained solutions are of applied interest, since they can be used for a better convergence 
of the original equation solution (§],|5|), which are true for arbitrary dimensions of the 
coupling hole. 
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7 Appendix 

If fi(z) can be representable as series (|13|), then 

Li(x) = / zJ 1 (xz)f 1 (z)dz = 2 V w n G n (x). 

where 



(A.l) 



Ji(xz) J (z)dz + A 2 



00 2 2 Ji(x2;) J (z)dz 



o A 2 - z 2 

00 J 1 (xz)J (z)dz 
X 2 n -z 2 



-GW(x) + \ 2 n GW(x). (A.2) 



It is easy to demonstrate (see [f|), that the integral G^ 1 ) at x > 1 is equal G^ = 1/x. 
Let's consider the integral G$p(x) at x > 1: 



Ji(xz)J (z)dz 



io A 2 - z 2 
Having made the substitution y = a;2, we obtain 

Jo A„ — cr|r 
where a = 1/x, < a < 1. Using the expansion 

OO 1 

J 1 (y) = 2yJ (y)J2 

m 



(A.3) 



(A.4) 



( |A.4| ) will take the form 

G( 2 )(x) = 2«E / 



2«E 



n,2 \2 _ \2 
m =l " L 



a 



00 yJo{y)Jo{ay)dy 
'Jo (A 2 - a V) (A*, - y 2 ) " 

yJ (y)J (ay)dy r°° yJ (y)J (ay)dy 



2a £ 



n (A 2 -«V) (A^-y 2 ) 
1 



n/ 2 \ 2 _ \2 
m=l " 'V 



ft 



2p(l) _ P(2) 
n m 



(A.5) 



Calculation of the integrals under consideration is based on such expansion (see||), 

Ti J Q (aj3z) 



4 Mz) 



[Uz)Y {(5z) - U(3z)Y {z)] 



E 



Jo (aP\ 8 ) Jo {P\ 8 ) 



{ J 2 (X s ) (z 2 - A 2 ) ' 



(A.6) 



where < a < 1, < (3 < 1. 
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Let's consider ( |A.6[ ) at z = X m /{3 

J (a(3X s )J ((3X s )(3 2 1 J («A m ) 



E 



or 

oo 



0£ 



\ Jl (X s ) (XI - f3 2 X 2 s ) 2 Ji (A m ) A m 
J (a/3A s ) J (/3A S ) /3A S J (aA r 



- {ttA s Jf (A,) /2} (XI - P 2 X 2 S ) Ji (A m ) A 

At (3 — > 0, we have 

p(2) 



yJo(y) J o( a y) d y _ Jo (aA m ) 

o (A m -£/ 2 ) Ji(A m )A m ' 
The value P« = can be obtained from ( |A.7|) at A, m — ► X n /a. Then 

G (2) M =2«V J ° (aAm) 

n 1 ' ^AK-* 2 vj* m Ji(x m y 



From the expansion 



J (az) _ i _ 2 ^ 2 ~ J (aA 



at 2 = A n /a we have 



Consequently, 



^l( z ) m =l A m Jo (A m ) ^ 2 - A m 

J (aA m ) 1 1 



E 



'j A m Jo (A m ) A 2 - a 2 X 2 m 2X1 ' 
G^(x) = — = -—. 



n n 



and 



n ~ XX 



Finally, 



POO J-^ 

Ia(x) = / zJ\(xz)f\(z)dz = 2 w n G 
Jo r - ^ 

which is what needed to be proven. 



,. x x) = 0, x > 1, 

n=0 
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